Third moments of conserved charges as probes of QCD phase structure 
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The third moments of conserved charges, the baryon and electric charge numbers, and energy, as 
well as their mixed moments, carry more information on the state around the QCD phase boundary 
than previously proposed fluctuation observables and higher order moments. In particular, their 
signs give plenty of information on the location of the state created in relativistic heavy ion collisions 
in the temperature and baryon chemical potential plane. We demonstrate this with an effective 
model. 



PACS numbers: 12.38.Mh, 25.75.Nq, 24.60.Ky 

Quantum chromodynamics (QCD) is believed to have 
a rich phase structure in the temperature (T) and baryon 
chemical potential (/ib) plane. Lattice QCD calculations 
indicate that the chiral and deconfinement phase tran- 
sitions are a smooth crossover on the temperature axis 
[l|, while various models predict that the phase transi- 
tion becomes of first order at high density The ex- 
istence of the QCD critical point is thus expected. To 
map these components of the phase diagram on the T- 
fiB plane is one of the most challenging and stimulating 
subjects which may be achieved by relativistic heavy ion 
collisions [1]. 

Various observables have been proposed for this pur- 
pose [3, H, S, 0]- Most scenarios suggested so far are 
concerned with fluctuations, such as those of conserved 
charges, momentum distributions, slope parameters, and 
so forth. For example, fluctuations of conserved charges 
behave differently between the hadronic and quark-gluon 
plasma phases, and may be used as an indicator of the 
realization of the phase transition [1, |^ . The singularity 
at the critical point, at which the transition is of second 
order, may also cause enhancements of fluctuations if fire- 
balls created by heavy ion collisions pass near the critical 
point during the time evolution 0] . Because of finite 
size effects and critical slowing down, however, such sin- 
gularities are blurred and its experimental confirmation 
may not be possible [1,0]. In fact, so far no clear evidence 
for the critical point has been detected in event-by-event 
analyses jSj]. Approaches to use higher order moments 
for this purpose have been also suggested recently 
and experimental attempts to measure those higher or- 
der moments were reported, for example, in Ref. llj . 



Almost all previous studies, however, focus on the abso- 
lute value, especially the enhancement, of each observable 
around the phase boundary. 

In the present Letter, we propose to employ signs of 
third moments of conserved charges around the averages, 
which we call, for simplicity, the third moments in the fol- 
lowing, to infer the states created by heavy ion collisions. 
In particular, we consider third moments of conserved 
quantities, the net baryon and electric charge numbers. 
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where Nc with c = B, Q represent the net baryon and 
electric charge numbers in a subvolume V, respectively, 
E denotes the total energy in V, 6Nc = Nc — {Nc), and 
6E = E — (E) . We also make use of the mixed moments 
defined as follows: 
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To understand the behaviors of these moments around 
the QCD phase boundary, we first notice that the mo- 
ments Eqs. (IlJ and ([2]) are related to derivatives of the 
thermodynamic potential per unit volume, lu, up to third 
order with respect to the corresponding chemical poten- 
tials and T. The simplest example is 7773 (BBB), which is 
given by 
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where the baryon number susceptibility, xb, is defined as 
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The baryon number susceptibility xb diverges at the 
critical point and has a peak structure around there 
[1, [13, [Hi- Since 7773 (BBB) is given by the ^b deriva- 
tive of Xb as in Eq. ([3]), the existence of the peak in xb 
means that 7r73(BBB) changes its sign there. Although 
the precise size and shape of the critical region are not 
known, various models predict that the peak structure of 
Xb well survives far along the crossover line [1, [H, 14 1 
(See, Fig. [T] as a demonstration of this feature in a sim- 
ple effective model; the details will be explained later). 
This means that the near (hadron) and far (quark-gluon) 
sides of the QCD phase boundary can be distinguished 
by the sign of 7773 (BBB) over a rather wide range around 
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the critical point. It is this feature that third moments 
carries more information than fluctuations (second mo- 
ments); fluctuations are, by deflnition, positive definite 
and cannot differentiate the near side from the far side 
as decisively as the third moments. We note that odd 
power moments around the averages in general do not 
vanish except the first order one. As we shall see later, 
all third moments presented in Eqs. ([T]) and ^ can be 
expressed in terms of derivatives of corresponding sus- 
ceptibilities which diverge at the QCD critical point, and 
hence change their signs there. 

The third moments can be measured in heavy ion col- 
lisions by the event-by-event analysis similarly to fluctu- 
ations, provided that Nc and/or E in a, given rapidity 
range, Ay, in flreballs created by collisions is determined 
in each event. The measurement of Nb is difficult be- 
cause of the difficulty in identifying neutrons. On the 
other hand, Nq and E can be measured with the exist- 
ing experimental techniques. Four out of the seven third 
moments in Eqs. H]) and ^ composed of Nq and E thus 
can be determined experimentally. 

All quantities we are considering here, A^b.q and E, are 
conserved charges and the variation of their local densi- 
ties requires diffusion. In Ref. Q , it was shown that the 
effect of diffusion is small enough for the fluctuations of 
the baryon and electric charges if the rapidity range is 
taken to be Ay > 1 [3] ■ In the estimate, the one dimen- 
sional Bjorken expansion and straight particle trajecto- 
ries were assumed. If the contraction of hadron phase due 
to the transverse expansion and the short mean free paths 
are taken into account, the above estimate will be more 
relaxed. This conclusion is not altered even if we take 
the effects of global charge conservation and resonance 
decays into account 0, 0, On the other hand, the 
experimentally measured charge fluctuations at RHIC ^] 
are close to those in the hadron phase than the free quark- 
gluon one. We, however, remark that the values of charge 
fluctuations in the quark-gluon phase can be similar to 
that in the hadron phase if the quark and gluons are 
strongly coupled 17[ ■ The experimental results therefore 
do not necessarily contradict a realization of the quark- 
gluon phase. 

Once the negativeness of third moments is established 
experimentally, it is direct evidence of two facts: (1) the 
existence of a peak structure of corresponding suscep- 
tibility in the phase diagram of QCD, and (2) the re- 
alization of hot matter beyond the peak, i.e. the quark- 
gluon plasma, in heavy ion collisions. We emphasize that 
this statement using the signs of third moments does not 
depend on any specific models. The experimental mea- 
surements of signs of moments also have an advantage 
compared to their absolute values: it is usually essential 
to normalize experimentally obtained values by exten- 
sive observables, such as the total charged particle num- 
ber A^ch, in order to compare the experimental results 
with theoretical predictions 0, 0] ■ In the measurement 



of signs, however, normalization is not necessary. In the 
measurement of absolute values, one has to be aware of 
the effect of global charge conservation when Ay is large 
[3]. It is, however, expected that the effect does not 
change the signs of the third moments. It is these fea- 
tures that our proposal is less subject to experimental 
and theoretical ambiguities and more robust than previ- 
ously proposed ones. 

Let us now consider the behavior of third moments 
other than m3(BBB) around the critical point. First, 
the third moment of the net electric charge 7713 (QQQ) is 
calculated to be 
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where /xq represents the chemical potential associated 
with Nq, i.e. d/d^iQ = {2/3)0/ d^iu - [1/3)8/ = 
((9/9/iB + d/dfii) /2, and the isospin chemical potential is 
defined as /zi = (/iu — /^d)/2 with /iu^d being the chemical 
potentials of the up and down quarks, respectively. In 
relativistic heavy ion collisions, the effect of isospin sym- 
metry breaking is small. Assuming the isospin symmetry, 
the second and last terms in the most right hand side of 
Eq. ([5]) vanish and one obtains 
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with the isospin susceptibility xi = —d^uj/dfj^f. Under 
the isospin symmetry, xi does not diverge at the critical 
point because the critical fluctuation does not couple to 
the isospin density [3] ■ The critical behavior of the term 
in the parenthesis in Eq. ([6]) in the vicinity of the critical 
point is thus solely governed by xb- Since TO3(QQQ) is a 
/iB derivative of this term, a similar behavior as 777.3 (BBB) 
is expected. 

Next, it can be shown that mixed moments including 
a single E are concisely given by 
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with c = B, Q, where xq = -d'^^/dfJ-Q = (xb + Xi)/4 
is the electric charge susceptibility. The T derivative in 
Eq. ([7|) is taken along the radial direction from the ori- 
gin with flxed fi = hb/T, i.e. d/dT\f, = d/dT\^^ + 
{fiB/T)d/dnB\T- Since Txc diverges at the critical point, 
Eq. ([7]) again leads to a similar behavior of 7773 (ccE) as 
the above-mentioned moments. 

To argue the behaviors of remaining third moments 
including two or three E's, it is convenient to flrst de- 
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moments are then given by 
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Since C/; is the second derivative of cu along the radial 
direction, it diverges at the critical point which belongs 
to the same universality class as that of the 3D Ising 
model. Therefore, m3(EEE), ma (BEE), and m3(QEE), 
all change their signs at the critical point. 

While the above arguments, based on the divergence 
of second derivative of oj, guarantee the appearance of 
the region with negative third moments in the vicinity of 
the critical point, they do not tell us anything about the 
size of these regions in the T-/xb plane. In fact, all third 
moments considered here become positive at sufficiently 
high T and /ib > where the system approaches a free 
quark and gluon system. The regions are thus limited 
more or less near the critical point. 

The information about the behavior of the third mo- 
ments at small /zb can be extracted from the numerical 
results in lattice QCD. For example, with the Taylor ex- 
pansion method the thermodynamic potential is calcu- 
lated to be w = -c2{T)nl - C4(T)^| - ce{T)n% , 

and one can read off the behavior of m3(BBB) at small 

fj.B asm3(BBB) = 24[c4(r)/iB + 5c6(r)/x| H ]. Lattice 

simulations indicate that Ci{T) is positive definite, while 
cg^(T) becomes negative in the high temperature phase 
[18|. From this result one sees that m3(BBB) is positive 
for small ^b, while the negative ce{T) suggests that the 
sign of m3(BBB) eventually changes at sufficiently large 
fiB- Other moments for small hb can also be evaluated 
in the Taylor expansion method by expanding w with 
respect to T and /xq. If the contour lines of vanishing 
third moments are close enough to the T-axis, the lattice 
simulations may be able to determine these lines. Since 
the region with a negative third moment should depend 
on the channel, combined information of signs of different 
third moments, and the comparison of the third moments 
obtained by experiments and lattice simulations, will pro- 
vide a deep understanding about the state of the system 
in the early stage of relativistic heavy ion collisions and 
the QCD phase diagram. 

The range of hb/T where lattice simulations are suc- 
cessfully applied, however, is limited to small ^b/T with 
the present algorithms. In particular, thermodynamics 
around the critical point cannot be analyzed with the 
Taylor expansion method. In order to evaluate the qual- 
itative behavior of the third moments in such a region, 
one has to resort to effective models of QCD. To make 
such an estimate, here we employ the two-flavor Nambu- 
Jona-Lasinio model [3, [ioj with the standard interac- 
tion £int = G{{'4)^)'^ + (V'*75'riV')^}j where tp denotes the 
quark field. For the model parameters, we take the val- 
ues determined in Ref. uM\ G = 5.5GeV~^, the current 




FIG. 1: (color online). T and /ib dependence of the baryon 
number susceptibility xb multiplied by T in the Nambu-Jona- 
Lasinio model. The bold line on the bottom surface shows the 
first order phase transition line and the point at the end is 
the critical point. 
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FIG. 2: (color online). Regions where third moments take 
negative values in the T-^b plane. The regions are inside the 
boundaries given by the lines. 



quark mass m = 5.5MeV, and the three-momentum cut- 
off A = 631MeV. For the isospin symmetric matter, this 
model gives a first order phase transition at large /zb, as 
shown on the bottom surface of Fig. [T] by the bold line. 
The critical point is at (T,^b) ^ (48,980)MeV. 

In Fig. [11 we also show the T and ^b dependence of 
TxB calculated in the mean-field approximation. One 
observes that xb diverges at the critical point, and the 
peak structure well survives along the crossover line up to 
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higher temperatures . The region where each moment 
becomes negative in the T-/LtB plane is shown in Fig. [2l 
One sees that ah the moments become negative on the 
far side of the critical point as it should be, whereas the 
extent of the region depends on the channel. The figure 
shows that areas with 7713 (BBB) < and TO3(BBE) < 
extend to much lower fiB and much higher T than the 
critical point. This suggests that even if the critical point 
is located at high /xb the negative third moments can be 
observed by heavy ion collision experiments. The figure 
also shows that the areas have considerable thicknesses 
along the radial direction. Since the system stays near 
the phase transition line considerably long regardless of 
the order of the phase transition, first order or crossover, 
once the state on the far side is created, negative third 
moments are very likely to be formed and observed. The 
wide regions of negative moments also indicate that they 
are hardly affected by critical slowing down and finite 
volume effects during the dynamical evolution of fireballs. 

Figure[2]also shows that areas with negative m3(EEE), 
m3(QEE) and 7713 (BEE) are much larger than those of 
the other moments in the T-fiB plane; although not 
shown in the figure, these areas extend even to the T- 
axis. The behaviors of 7713 (EEE) and m3(cEE) near the 
T-axis can be checked directly by the lattice simulations. 
If the range of T satisfying 7713 (EEE) < is sufficiently 
wide at /ie = 0, it is possible that the negative third 
moments are measured even at the RHIC and LHC en- 
ergies. Whether the negative moments survive or not 
in this case depends on the diffusion time of the energy 
density, in other words the heat conductivity. One can 
thus use the signs of 7773 (EEE) and 70,3 (cEE) to estimate 
the diffusion time of the charges and energy. The third 
moments 777,3 (QQQ) and 7773 (QQE), on the other hand, 
become negative only in small regions near the critical 
point. These behaviors come from the large contribution 
of XI in Eq. dSD. 

It should be, however, remembered that the results in 
Figs. [Hand m are obtained in an effective model. In par- 
ticular, the model employed here gives the critical point 
at relatively low T and high /iB [2| • If the critical point 
is at much lower /iB , the areas with negative moments in 
Fig. [5] should also move toward lower and higher T. 

In this Letter, we have pointed out that the third mo- 
ments of conserved charges, the net baryon and electric 
charge numbers and the energy, carry more information 
on the state around the QCD phase boundary than usual 
fluctuation observables. They change signs at the phase 
boundary corresponding to the existence of the peaks of 
susceptibilities. If the negative third moments grow at 
early stage of the time evolution of fireball created in the 
collisions and if the diffusion of charges is slow enough, 
then the negative third moments will be measured exper- 
imentally through event-by-event analyses. Once such 
signals are measured, they serve direct evidence that the 
peak structure of corresponding susceptibility exists in 



the phase diagram of QCD, and that the matter on the 
far side of the phase transition, i.e. the quark-gluon 
plasma, is created. The combination of the third mo- 
ments of different channels, and their comparison with 
the numerical results in lattice QCD will bring various 
information on the phase structure and initial states cre- 
ated in heavy ion collisions at different energies. 
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